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Spatially extended defects such as the F center give rise to an effective defect charge density 
which may produce important polarizations in the crystal. The electric field in the crystal depends 
upon these induced dipoles. Lattice summations for the contribution to the electric field which arises 
from ionic shells centered about the defect are evaluated for the NaCl and CaF 2 lattice structures. 
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1. Introduction 

We shall evaluate in this paper the ionic and elec- 
tronic polarizations which defects induce in ionic 
crystals. Such polarizations play a major role in the 
energy of formation, the mobility, the self energy, and 
the electronic structure of defects. 

We shall consider an ionic crystal which contains 
one molecule M n X„ per volume v c of the crystal. 

The volume v c is the volume of a unit cell and is given 
by 



and at 



v c = ai ' (a 2 X a 3 ), 



(1) 



where ai, a 2 , and a 3 are the primitive translational 
vectors. For NaCl and CaF 2 structures, the primi- 
tive translational vectors are a t =(a/2)(0, 1, 1), 
a 2 =(a/2)(l, 0, 1) and a.,= (a/2)(l, 1, 0), and the 
volume of the unit cell is v r = a ; 74. The lattice constant 
a is the cation-cation distance. The nearest neighbor 
distance (cation-anion) is r*o= («/2) for NaCl structures 

and r {) = ( v3/4)a for CaF 2 structures. The cations are at 

r + = (a/2)[(/ 2 + k )x + (/:, + /, )y+ (/, + 1 2 )«] . 
The anions for NaCl structures are at 

r- = r++(a/4)(i + f+z), 
and the anions for CaF 2 structures are at 



r, =r^ 



(a/4)(£ + y+i), 



r 2 - = r + -r- (a/4) (— x — y — z) . 



The integers // have the values 0, ±1, ±2, etc. The 
above discussion places a cation at the origin. We 
shall treat also the case for which an anion is at the 
origin. For this latter case we have the reference 
vectors r a =(a/2)(0, 0, 1) for NaCl structures and 
r (l = (a/4)(— 1, —1,-1) for CaF 2 structures. 

We view the anions X and the cations M as polariz- 
able point charges occupying the lattice sites. The 
anion ionicity is Z a and the cation ionicity is Z c . The 
unit cell is electrically neutral, n + Z c J rn-Z (l = 0. We 
also associate an ionic polarizability ofi a or a** and an 
electronic polarizability a e a or a e c with each ion in the 
lattice. These polarizabilities determine the response 
of the crystal to weak static and high frequency elec- 
tric fields. We represent an ion which is located at the 
lattice site r v and which experiences a static electric 
field E(r„) as a point charge Z (l e orZ r e located at r p and 
upon which we superpose a point dipole moment 



li(r„)=(a e +a d )E(x v ), 



(2) 



When the period of the electric field is much shorter 
than the characteristic period for ionic motion but also 
longer than the orbital period for the electrons about 
the ion, the point dipole moment is 



IJL(r v )= a e E(r„). 



(3) 



We assume that the response is linear in the electric 
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field. In general, the polarizability is a tensor. However, 
we choose here for convenience that it be diagonal 
and isotropic. This is a very reasonable assumption for 
cubic crystals. When we view an ionic crystal as a 
dielectric continuum, the static dielectric constant 
€o and the high frequency (optical) dielectric constant 
e x give us information about the ionic and electronic 
dipoles which are induced by weak electric fields. 

For our purposes here, we may represent the defect 
by a spatially dependent charge density p(r) [I]. 1 
For example, the F center in an ionic crystal consists 
of an electron (the F electron) localized about an anion 
vacancy. We choose the anion vacancy to be at the 
lattice site r = 0. We characterize quantum mechan- 
ically the F electron by its stationary wave function 
i//»(r) with energy eigenvalue e n where n represents all 
the quantum numbers of the given state. The F elec- 
tron gives rise to the spatial charge density 

p r (r)--ei//*(r)i//„(r). 

Also, we view the anion vacancy as the addition of a 
positive point charge Z v e = — Z (l e at the site r = 
p v (r) =Z v ed 3 (r) . The total charge density pd(r) of the 
F center at r = then becomes. 



p d (r)=e[Z r 8 3 (r)-^*(rW ll (r)]. 



(4) 



The total electric field which a polarizable ion 
experiences is the sum of the electric field due to the 
defect charge density p<f(r) and the electric field due 
to the electronic dipoles and the ionic dipoles which 
the defect charge density induces on all the other ions 
of the crystal. A. D. Franklin and D. J. Sparks [2] have 
computed the electric fields produced by a point charge 
p p (r) = Z eff eS 3 (r) in ionic crystals. We will compute 
in this paper the electric fields produced by extended 
charge distributions in the NaCl and CaF 2 lattice 
structures. 

Our intended research on the electronic structure 
of the F center motivates us to carry out the present 
computations. The F electron and the lattice, par- 
ticularly the nearest neighbors to the defect, form a 
very strongly interacting system which we must treat 
self-consistently. The polarizations which accompany 
the F center are most important when we use the 
polarizable ion theory for the F center [1]. 

2. Formulation 

In order to make the present computations feasible, 
we shall consider only the spherically symmetric part 
p s (r) of the charge density pd(r) ; i.e., we average p d (r) 
over the unit sphere. 



p.,(r) = (477)" 



r 27r r u 

d<t> 

Jo Jo 



sin 4>d4>p(r) 



(5) 



"or convenience, we will suppress 



D 



tin 



where r= |r| 
subscript s. 

The electrostatic potential </>( r i) at the point ri due 
to the spherically symmetric part p(r) of the defect 
charge density prf(r) is given by the expression, 



<t>(n) = 



p(r)r*dr 

Jo r pjrjr^dr 

n J n r 



(6) 



We do not allow the dielectric response of the system 
to enter in eq (6). The negative gradient of the elec- 
trostatic potential then prescribes the electric field 
Ee/(ri) at r t due to the charge density p(r); 



E f /(ri)=£ rf (ri)r, = - 



d<Mr,) 



Jo 



p{r)j^dr 



A 



(7) 



where ri = ri/|ri| is the unit vector in the direction n. 
We also define the average amount of charge Q(r) con- 
tained within a sphere of radius r centered about 
r=0, 



<?(r) 



Jo 



{rMdr,. 



Combining eqs (7) and (8) gives us 
E d (r) = r-*Q(r). 



(8) 



(9) 



The electric field E f /(r,,) at rv due to the defect 
charge density induces a dipole moment fl(r v ) on the 
ion at site r„. The electric field at the site iv, due to the 
dipole moment fi(r v ) at the site v v is 



-■^dipl'Vi? *V 



== 3(r, 1 -r r )[/x(r,) 

Ir^-rv 



(r^-iy)] 



liVi — ivl 



(10) 



i?tot(r„)=Ed(r„)+ ^ E di P (iv; iv,), 



(11) 



1 Figures in brackets indicate the literature references at the end of this paper. 



where the prime denotes a sum over all lattices sites 
except the site y Vx = y v . Multiplying the left and right 
hand sides of eq (11) by the appropriate polarizability 
for the ion at r„ converts eq (11) to a system of 3/V equa- 
tions having 3./V unknown variables — the three com- 
ponents of the dipole moments §Ji(r v ) at each of the /V 
sites in the lattice. Because /V has the same order of 
magnitude as Avogadro's number, we cannot solve 
the system of equations given by eq (11). Hence, we 
must make further approximations. The approxima- 



The total electric field which the polarizable ion at 
r„ experiences is the sum of the electric field due to the 
defect and the electric field due to the dipoles induced 
on all the other ions except the ion at i\, 
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tions which we will employ to overcome the above 
formidable task are known as the Mott-Littleton 
procedure [3], 

We divide the crystal into two regions. Region I 
contains all the ions on shells centered about the defect 
and having a radius less than or equal to r s . The radius 
r,s is the radius of the sth shell of ions centered about 
the defect. Region II contains all the remaining ions. 
The Mott-Littleton procedure solves the system of 
equations given by eq (1 1) for only region I and employs 
a continuum model to approximate the solutions for 
the dipole moments of region II. 

We first outline the prescription for the dipole mo- 
ments in region II. We view region II as a dielectric 
continuum. The electric field E and the displacement 
vector D are related in region II by the constitutive 
equation for a dielectric continuum, 



E„ u ,(r 1 )= 2 



v QJj. 


)M n {v) 




f ^|r,-iv| 3 




V* r « 






"3(ri — r„){iv -{r, 


-r„)}-r„" 




[ |r,-iv| 2 





D(r)=eE(r)=E(r)+47rP(r), 



(12) 



where e is the appropriate dielectric constant and the 
polarization P(r) is the dipole moment per unit volume, 

P(r)=[(e-l)/47r]E(r). (13) 

The dipole moment per unit cell at n is, 

P c (r,) = v c P(r) = Vc[(e- 1 )/4rr|E(r) (14) 

The Mott-Littleton presecription asserts that all dipole 
moments point in a radial direction from the defect 
and divides the dipole moment per unit cell among 
the n + + Ji- ions contained in the unit cell in proportion 
to their individual polarizabilities; 



(17) 



The summation is over all the sites in region II. The 
radial part of E out (r,), £ out (n) = E out (n) • (pi/|pi|), is 
given by, 

E out (ri)= Y — rr 

*f rir/,|r, — 1\,|° 



X [2(r, • r,)r?-3^- (r„ • r,) 2 + 2rf,(r, • r,)]. (18) 

When we express the summation which appears in 
eq (18) in terms of a summation over the shells, we may 
evaluate more easily the radial part of E out (ri). Each 
shell of ions centered about the defect has ions of the 
same type and each ion of shell s exhibits the same 
polarizability M u (v) =M n (s). Hence, we may write 
an alternate expression, which is computationally 
more practical, for the radial field E i)[lt (r\); namely, 

Eout(r,) = r;r 8 £ QMM n (*')DU', r,). (19) 



The dipole coefficients for each shell s, which is out- 
side region I, D (s; Pi) is a summation over the n (s) ions 
on the .sth shell; 



n(s) 

D(5;n) = 2 



n 



nr;|ri 



fin(r v )= 



a v 



{n+OLc + n-a (l 



v c 



(e-l)E„(p,.) 
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€1 



fr (15) X [2(r, w • ri)r?-3/trf- (p. • r,) 2 + 2r?„(r„ 



i-i 



)]. (20) 



When treating the F center problem one must distin- 
guish between the electronic polarizations and the 
ionic polarizations. The electronic polarizations re- 
spond to the rapid changes in pd(r) which occur when- 
ever the F electron undergoes a transition from one 
state to another state. However, the ionic polarizations 
may not respond to such rapid changes inp f f(r). Hence, 
when ei = e x , then e must be 6 X . But when €i = €o, 
e = €oc when we compute only the electronic polariza- 
tions and € = €o when we compute both the electronic 
and ionic polarizations. When both ionic and elec- 
tronic polarizations occur a,,= a 
electronic polarizations occur a,. 
(15) we define the polarizability of an ion in region II 
by the quantity. 



-a^ and when only 
a e u . Referring to eq 



M u (v) = 



a r vAe - 



^TT€\{n + a v + n-(x (i ) 



(16) 



The electric field at a site r, in region I due to all 
the dipoles [eq (15)] in region II is 



where the site r m is the position of the with ion of shells 
and |r,„| = r,.. 

We may show from cubic symmetry arguments and 
from the assumption of radial dipoles that the dipole 
coefficient D(s, ri) evaluated for the lattice site ri on 
shell t in region I is the same for all the other lattice 
sites r m on the same shell t in region I. 

The system of equations for the dipole moments, 
eq (11), is tractable only when region I contains a small 
number of shells; e.g., u shells. Incorporating eqs 
(19) and (20) into eq (11) and assuming that the dipole 
moments in region I also are radially directed, we may 
reduce the system of equations given by eq (11) to u 
equations containing the u dipole moments on the ions 
in region I. The magnitude of each dipole moment at 
the lattice sites on a given shell will be the same. 

3. Results 

The experience which we have gained from our 
study of the F center in the alkaline earth oxide CaO [4] 
indicates that reasonable agreement with experiment 
obtains when we use the first order Mott-Littleton pro- 
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cedure to evaluate the polarizations associated with 
the F center. We shall present therefore the first order 
Mott-Littleton prescription in order to demonstrate a 
specific example of how one might use the preceding 
discussion. 

The first order Mott-Littleton approximation states 
that region I contains only the first shell and that 
region II contains all the remaining shells. The first 
shell about an anion defect contains n c cations at a 
distance r\ — r () from the defect center. Again, we as- 
sume that the induced dipole moment on an ion in the 
first shell points in a radial direction and has a magni- 
tude /jl\. The radial electric field at the site ri due to the 
dipole moments /Xir™ on the (n c — 1) other ions on the 
first shell is 

ri ' 

where the same shell dipole coefficient C\ is a constant 
for a given lattice structure. The system of equations 
represented by eq (11) reduces in the first order Mott- 
Littleton approximation to one equation for the one 
unknown dipole moment jUi; namely, 



where E out (ri) is given by eq (19) in terms of the dipole 
coefficients D(s, ri). 

We list in tables 1 and 2 the dipole coefficients 
D{s, ri) for the NaCl and CaF 2 structures. The com- 
puter program used to evaluate the dipole coefficients 
may be easily modified to treat defects centered about 
cations and to treat any number of shells in region I. 
The results given in tables 1 and 2 are applicable 
directly to the F center problem. 



Em=- 



(21) 



Ml : 



'Q(ri) 

ri 



+£out(ri) + 



MiCi 



(22) 



Equation (22) then predicts that the dipole moment ix\t\ 
on the ions in the first shell about the defect is 



flih- 



L r\ 



1 — 



nOt\ 



r\ j 



fu 



(23) 



Table 1. The dipole coefficients D(s, ri) for NaCl structures 

The defect site is an anion site and the evaluation site ri is a cation site on the first nearest 
neighbor shell centered about the anion defect site. Region I contains only the first nearest 
neighbor shell. Region II contains all the remaining shells. The number of cations in the 
first shell is n c — 6 and the same shell dipole coefficient C\ for the (n c — 1)= 5 other ions on 
the first shell is C, = — 2.3713. The radius of shell 5 is r, in units of (a/4) and the dipole co- 
efficients D(s, ri) for the shells outside the first shell are given by eq (20). The lattice constant 
is a and the nearest neighbor distance is r n = (a/2). The dipole coefficients are expressed in 
the form Vx 10". 
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The author thanks A. D. Franklin for many informa- 
tive discussions. 



Table 2. The dipole coefficients D(s, ri) for CaF 2 structures 

The defect site is an anion site and the evaluation site ri is a cation site on the first nearest 
neighbor shell centered about the anion defect site. Region I contains only the first nearest 
neighbor shell. Region II contains all the remaining shells. The number of cations in the 
first shell is n r — i and the same shell dipole coefficient C\ for the (n r — 1)=3 other ions on 
the first shell is d = — 1.1482. The radius of shell s is r x in units of (a/4) and the dipole co- 
efficients D{s, r, ) for the shells outside the first shell are given by eq (20). The lattice constant 
is a and the nearest neighbor distance is r = (3 1/2 a/4). The same shell coefficient d given 
in ref. 2 is expressed in terms of the anion-anion distance and the above value is in terms of 
the anion-cation distance. The dipole coefficients are expressed in the form /VX 10". 
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2#<s, r,) 
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